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We present a calculation of the zeta function and of the functional determinant for a Laplace-type 
differential operator, corresponding to a scalar field in a higher dimensional de Sitter brane back- 
ground, which consists of a higher dimensional anti-de Sitter bulk spacetime bounded by a de Sitter 
section, representing a brane. Contrary to the existing examples, which all make use of conformal 
transformations, we evaluate the zeta function working directly with the higher dimensional wave 
operator. We also consider a generic mass term and coupling to curvature, generalizing previous 
results. The massless, conformally coupled case is obtained as a limit of the general result and 
compared with known calculations. In the limit of large anti-de Sitter radius, the zeta determinant 
' for the ball is recovered in perfect agreement with known expressions, providing an interesting check 

' of our result and an alternative way of obtaining the ball determinant. 
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' The idea of living on a membrane embedded in a higher-dimensional spacetime has attracted enormous attention 
, over the past few years, the main motivation being the fact that the localization of particles on branes provides an 

alternative to the standard picture of Kaluza-Klein compactification P, |3| . 
' A popular example is the Randall-Sundrum model 3], which considers a five-dimensional slice of Anti-de Sitter 
(AdS) spacetime with the extra dimension compactified to an orbifold having two 3-branes of opposite tension located 
at its fixed points. This results in a non-factorizable spacetime, which, as Randall and Sundrum have shown, has 
the important consequence that the effective four-dimensional scale on the negative-tension brane turns out to be 
O ^' exponentially suppressed relative to the higher-dimensional scale. This was originally proposed as a solution to the 
QJ , hierarchy problem by explaining the very small ratio (large hierarchy) of some 17 orders of magnitude between the 
i-C ' electroweak scale and the Planck scale observed on our brane, identified as that with negative tension. In this model, 
the hierarchy arises as a geometrical effect, with gravity being strongly localized around the positive-tension brane. 

Pushing this idea further, Randall and Sundrum |^ showed that the negative-tension brane may actually be absent 
and the extra dimension possibly infinite in size. In that case, we live on the positive-tension brane with gravity 
^ . confined around it, recovering its standard Einstein form from an effective four-dimensional point of view. This result 
does not solve the hierarchy problem, but has stimulated the construction of many new cosmological models. Here 
we concentrate our attention on the de Sitter brane model, relevant in the construction of brane world models of 
infiation. 

The problem of studying quantum effects in such scenarios is, then, naturally posed. It is immediately clear that 
when a quantum field is considered on such background spacetimes, quantum effects may play a significant role, the 
simple contribution that they give to the brane and bulk cosmological constants being an evident example. 

A number of people, inspired by previous work in Kaluza-Klein theories [^|^, have investigated the possibility that 
quantum effects from bulk fields could play a role in the consistency of such models by providing a sensible mechanism 
of stabilization in the Randall-Sundrum two-brane model and in some higher-dimensional generalizations. Various 
authors have dealt with the problem of quantum effects in brane models and some references are |^ 1^ lol lTollTllfl^ 
IT^Ilslfl^ . The lowest order quantum corrections arising from bulk fields on the Randall-Sundrum background have 
been calculated in a variety of ways, and then extended to some general classes of higher-dimensional spacetimes in 
Refs. Ref. |23| considers the case of scalar and gauge fields in the Randall Sundrum model and interprets 

the result in terms of the AdS/CFT correspondence, showing explicitly that quantum effects could provide a sensible 
stabilization mechanism. 
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It is worth mentioning that in most of the previous work, the geometry of the branes was assumed to be flat, which 
greatly simphfies the study of quantum effects. When the branes are curved, as happens, for example, in de Sitter or 
hyperbolic brane models, and where the Casimir energy could have some effect on the cosmological evolution of the 
brane, the situation becomes more complicated. 

Some work in this direction has been carried out in Refs. [2lll2^ I2M |2^ . Specifically, in Refs. 0,112, the vacuum 
energy for a massless conformally coupled scalar field in a brane world corresponding to de Sitter branes in an AdS 
bulk has been evaluated. This calculation, which is technically the simplest, is carried out by working in a conformally 
related spacetime, similar in form to the Einstein universe. Zeta function regularization is employed and the result 
shows that the vacuum energy is zero for the one brane configuration. These results have been extended to the case 
of conformally coupled Majorana spinors in Ref. j^J], still by making use of conformal transformations. The cases of 
a massless conformally coupled scalar field and of a massive minimally coupled field for de Sitter branes embedded 
in both de Sitter and AdS bulks have been considered in Ref. p3l | , where the vacuum energy is computed once a gain 
by using conformal and zeta regularization techniques and is found to be zero for the one brane case. Refs. [21I l25j| 
consider a somewhat related set-up and compute the effective action for scalar fields in an AdS bulk bounded by AdS 
branes, still by making use of conformal transformations. 

Apart from the relatively simple case of massless, conformally coupled field, the general case has not yet been 
studied. The present paper is devoted to providing a new derivation of the zeta function and of the functional 
determinant for a scalar field in a de Sitter brane model. The method is very general and applies, in principle, to 
a number of situations, where the methods based on conformal transformations or dimensional reduction cannot be 
applied at ease. 

We perform the calculation by directly working in the higher-dimensional spacetime and evaluate the zeta function 
for the original higher-dimensional wave operator. With respect to the previously studied cases, we have to deal with 
two main problems. The first is to compute a zeta function where the operator spectrum is not known, for which 
we adopt the technique developed in Refs. |2^ ll^, IH, |2^ [s^, |^ 113 • This general technique, elaborated in various 
forms, allows one to obtain the zeta function using only the knowledge of the basis functions. Such a method is 
applicable whenever an implicit equation satisfied by the eigenvalues is known, and has already been applied in the 
context of brane models. We stress, however, that in the case of flat branes the basis functions can be conveniently 
expressed in terms of Bessel functions, which greatly simplifles the calculation. The second main problem is that, 
for the class of spacetimes we consider, namely an AdS bulk bounded by a de Sitter brane, the basis functions are 
expressed as a combination of generalized Legendre functions, which are considerably less manageable. However, the 
method developed in Ref. [13 proves to be particularly useful and we closely follow their approach in our calculation. 

The structure of the paper is as follows. In the next section, we discuss the spacetime configuration of the de Sitter 
brane model and solve the higher-dimensional scalar wave equation on such a background spacetime. For the sake of 
generality, we consider a massive scalar field coupled to the higher-dimensional curvature. In the subsequent section, 
after having described the general technique of Ref. [13, we pass to the main task at hand, which is the evaluation of 
the zeta function for a de Sitter brane model. After having obtained the general result, we consider two limiting cases, 
which provide a non trivial check on the method used as well as the actual calculation. First we specify the result 
to the conformally coupled case and compare with that of Refs. ^21, 22, 23, 24]. Then we consider the limit of large 
AdS curvature radius which should reproduce a ball-like geometry. The last section is left for conclusions. Several 
technical results regarding the asymptotic expansion of the generalized Legendre functions, as well as the derivations 
of certain results used in the calculation, are reported in the appendices. 



II. SCALAR FIELDS ON DE SITTER BRANE BACKGROUNDS 



where 



The background configuration we consider is described by the action 

S = J d^xV^ y^{n - 2A5)) - J d^xV^ T , (2.1) 

A5 = -| (2.2) 

is a higher-dimensional, negative, bulk cosmological constant term, r is the brane tension and is proportional to 
the five dimensional gravitational constant. The corresponding solution of the setup, with the usual cosmological 
symmetries, has an AdS five-dimensional geometry; £ is the radius of the five dimensional AdS space. 
A convenient scaling of the time coordinate allows one to write the metric as 

ds^ = dr^ + [Hlf smh^ {r/l)[-dt^ + H~^e^"*dyil] , (2.3) 

where the coordinate r parameterizes the extra dimension and the time coordinate corresponds to the cosmic time 
parameter on the brane. For notational convenience, we define 

a(r) = iJ£sinh(r/£). (2.4) 
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A Z2-symmetric brane world can then be constructed in a standard way by taking two slices of the space and gluing 
them along the brane located at vq. In such a case, the junction conditions at the brane give the Friedmann equation 



As 
6 



(2.5) 



In the present case, the Hubble parameter, H, is constant, so that the brane geometry is de Sitter. The Hubble 
parameter is related to the brane position vq according to 



1 



£2 sinh^(ro/^) 



(2.6) 



On such a background, we consider a bulk scalar field and ignore its back reaction. It satisfies the Klein-Gordon 
equation 



(2.7) 



For the sake of generality, we consider the brane to be D- rather than four-dimensional, for the remainder of the 
present section. Furthermore, for later convenience, we transform to Riemannian signature. Hence, De is the {D + 1)- 
dimensional d'Alembertian on Riemannian space and TZ is the scalar curvature, given by 



n 



2D 



a{r) 



D{D-1) 



a{r) 



(2.8) 



where TZs = D{D — l)iJ2 is the scalar curvature of the de Sitter brane, which is a D— sphere of radius H ^. 

We are interested in finding the eigenmodes (pnj{x, r) and eigenvalues ^ of the above field operator, defined by 

(-□e + +e7^) ^n.j{x,r) = A^^- ^n^j{x,r) . (2.9) 
Let us assume that the modes are separable in the variables x and r: 

^n,j{x,r) = (j)j{x)fn,j{r), (2.10) 

where the spherical eigenfunctions satisfy 

-U^<i>,{x)^3{i + D-l)H^,{x), (2.11) 
where is the d'Alembertian on the de Sitter section and dj the degeneracy factor. 



(2.12) 



with i = 0,1,2,.. .. Using Eq. 12.1111 in Eq. (|2.9|) allows us to find the equation of motion for the radial eigenfunctions: 



[-a-'^idra'^dr) + a-^j {j + D - I) + + ^U] /„,,(r) = A^ , /„,,(r) , 
whose solution can be written in terms of toroidal Legendre functions 

fnAv) = {Hey^^-^/^inh'/'-^/'irj) [A„,,i^;;;^^_^/2(coshry)+i?„,,Q-^';;;^._^/2(cosh'7) 
where we have defined t] — r/£. The order and degree of the associated Legendre functions are 

D-1 



^3 = \:i 



and 



■Jn,3 - \l i'^^lj 



where 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



Regularity at the origin implies that Bn,j 
alized Legendre functions. 

Thus, our eigenfunctions take the form 



0, as follows from examining the small- argument behavior of the gener- 



(2.17) 
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with 0(77) = Hlsinhrj. 

As well as the solution in the bulk, we must also consider the boundary condition on the brane, which can be 
obtained by integrating l|2.13(l across the brane. In general, the Z2 symmetry allows us to choose either an untwisted 
field configuration, such that /(jy) = /(—??), corresponding to Robin boundary conditions 

d^fnjivo) = ^coth(77o)/„j(??o) , (2.18) 

or alternatively a twisted field configuration, such that f{rj) = — /(— corresponding to Dirichlet boundary conditions 

fnAvo) = 0. (2.19) 
In the remainder of the paper we focus our attention to Dirichlet boundary conditions. 



III. ZETA FUNCTION FOR DE SITTER BRANE MODELS 



A. General method 



The main scope of this paper is to compute the zeta function and the functional determinant for a bulk scalar field 
on a de Sitter brane background. For calculational simplicity, we take a field obeying Dirichlet boundary conditions, 
although the method can be applied to more general situations with few calculational modifications. To deal with the 
case of de Sitter brane models, we follow the approach of Refs. [2fill2^ IsoL Isills^ . where a calculational technique for 
(■-functions of differential operators on manifolds with boundaries with explicitly unknown spectra has been developed 
and applied (see Refs. 113, IH,!!^) to evaluate the one-loop contribution from the graviton, and matter fields, to the 
Hartle-Hawking wave function. Here we summarize their method. 

It is a well-known fact that the one-loop effective action can be written as 

rfi^^^TrlnA, (3.1) 

where A is a second-order differential operator, in our case given by 

A = -DE + m^ +^TZ . (3.2) 
r^^) can be expressed in terms of a generalized (-function 

C{s) = Y,X~\ (3.3) 

A 

with A being the eigenvalues of the operator A, which we assume to be positive definite. One has 

rW^-ic'(0)-ic(0)lnA.2 , (3.4) 

where /i is the renormalization scale. Thus, we see that the main problem is reduced to that of evaluating the 
(■-function and its derivative at s = 0. 

Usually, the computation of the (-function requires explicit knowledge of the spectrum. However, in many situations 
of interest, the eigenvalues are not explicitly known. To bypass this kind of problem, various authors have developed a 
calculational technique that allows one to evaluate the (-function and related quantities like functional determinants, 
Casimir energies and effective actions, when such information on the eigenvalues is lacking and the only knowledge of 
the spectrum is via an implicit equation [27L Isol IsiL l3^ . Generally, one considers the eigenvalue problem associated 
with the operator A: 

Aui,{x) = \ui,{x) , (3.5) 

where the parameter 1/ enumerates the independent solutions of (|3.5|l . A degeneracy factor d{i') is associated with 
each V. Imposing the relevant boundary conditions leads to an equation of the form 

F(A,i^) = 0, (3.6) 



where the function F{X, v) depends on the mode functions, the eigenvalues A, the index u and eventually other 
parameters inessential for the present discussion. 
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To evaluate the C-function it is not necessary to solve the previous equation, as is clear by making use of the residue 
theorem, which allows one to write 

C(s) = f ^dMlnF(c^,;.)da; , (3.7) 

where the contour C is chosen to enclose all the positive solutions of H3.6|l in the complex ijj— plane. 

For the explicit calculation, it is convenient to express this contour integral as an integral along the real line, 
which can be achieved by appropriately deforming the contour C. Typically, if the function F(uj,i>) satisfies certain 
properties, as in the case we will consider in the next section, a choice of a contour like the one in Fig.^ allows one 
to rewrite the integral H3.7|l as 

Cis) = '^ r z-^^Y.d{^.)\nF{z,v)dz. (3.8) 



The source of divergences in the above expression come from the large v behavior and the integration over z. We 
need to check that large values of s indeed regulate C,{s). Let us first presuppose that this is possible and, assuming 
s to be large enough, we swap the order of integration and summation and consider the asymptotic behavior of the 
integral 

z-'-^\nF{z,v)dz (3.9) 
dz 

for V —t oo and its integrand for z — > cx3 and z ^ 0. Now, it is possible to obtain a uniform asymptotic expansion 
$(i^, z/v) of lnF(z, v) such that v oo, while the ratio z/v is held fixed. Thus 

\nF{z,v) ^ <^{v,z/v) , (3.10) 

for large v. There are two important features of this expansion, aside from its uniformity: it has a power-law behavior 
of fixed order in v, as we shall illustrate for our case, and it is valid in the full range of the ratio < z/v < oo. If 
the integral 1)3.911 is finite for a large enough value of s, which is true in our case, then uniformity, together with these 
two properties ensures the convergence of the sum over (|3.9ll . for sufficiently large s. This can be proved under quite 
general assumptions, but here we content ourselves with showing that this is true for the case under study. Finally, 
a simple rescaling allows us to write 

C(^) = ^^d(^)^- r z-^^<^{z,v)dz . (3.11) 

TT ^ Jo dz 

This verifies that it is possible to regularize the divergences by a suitable choice of s. We now interchange the 
integration and summation back to their original order and write 



^, , smTTS 
C(s) = 



r°° d 

/ z-''—I{z,s)dz, (3.12) 
Jo dz 



where I{z, s) is given by 



/(z,s) = ^d(i/)i/-'*-^lni^(z,i/) . (3.13) 



In order to further proceed with the evaluation of the ^-function, we must analytically continue to s = 0. We 
expand the sum H3.13|l around small values of s, which generally develops a pole: 

I{z,s) Li+/R(^)+o(s) . (3.14) 

s 

It is now possible to use the following lemma along with the properties of the asymptotic expansion to compute 
C(s). Consider a function f{x) which is analytic at a; = e, for some small e, and has the following general asymptotic 
behavior, when a; — > oo, 

Pk<N 

fi^)^ ifk + ~fk\r,x)x'"' + \f]io^h-ix+[fUs + Oix-') , pk>0, (3.15) 

k=l 
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FIG. 1: The original contour (enclosing the positive real zeros) and the deformed contour lying along the imaginary axis. The 
deformed contour C' is closed at infinity by a semi-circle of infinite radius. 

where the subscripts log and reg refer to the solely logarithmic and regular (non-singular) parts oi J(x) in the large 
X limit. Then, there exists the analytic continuation of the integral 

where [f\f = [f]rc^ - /(e), for example see Ref. p7| . 

On the basis of the uniform asymptotic expansion of the eigenfunctions, it is possible to prove that /p°'°(0) and 
/^ (z) behave as (|3.15|) and we also assume that [/P°''^]iog = 0, which we shall show to be true in the case of interest 
to us. It is now a simple matter to apply the result H3.16|l to /p°'''(z) and /^(z) to get 

C(s) = niog + [/P°'°]o°° + s|[^%°°-y^ dzlnz2^_U|+o(,2) ^ ,^0. (3.17) 
This equation can then be used to get the value of the ^-function and its derivative at s = 0. 

B. Explicit evaluation of the ^-function 

We now pass to the explicit evaluation of the ^-function for the scalar field on the de Sitter background, described 
in section (QfJ. We have seen that, for Dirichlet boundary conditions, the eigenvalues are given by the solution of the 
implicit equation 

F{tOnS. ^i) = 2'^r(l + Vj)P72^_^i^[cos\iJ]^) = . (3.18) 

Note that we have multiplied the Legendre function by the factor -I- Vj) as this does not change the solution 

of the eigenvalue equation H2.19|l and, on the other hand, produces some simplifications in the calculations at a later 
stage. It is also clear that any factor independent of z does not affect the contour integral 1)3. 7|l . 
Thus the ^-function can be expressed by the double sum 

c(s) = = E K.. + ^')"'' (3.19) 

n,j n,j 

with LOn,j defined in (|2.15|l being the solutions to (|3.18|) . As described in the previous subsection, we can use the 
residue theorem to express the ^-function as a contour integral, in terms of the complex parameter co 

p2s r J J 
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where the contour C is chosen to enclose the real positive zeros of F{uj, Vj). 

By appropriately deforming the contour (see Fig. ^) and by performing some formal manipulations, we arrive at 

a-^) = ^l ^d^n^^s). (3.21) 

where 

/(z2,5)^^d,-^ ln[2''^T(l + ^.,) p-^_,_^_^(coshr;o)] . (3.22) 

7 V J 
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J 



As in the previous subsection, we split the contributions to the function into one regular plus one polar piece: 

C(s) = Cfl(-?)+Cp(s) , (3.23) 

where 

Us) = !^ r d ^i^(^2^ ^) (3 24) 

TT Jq az 



and 



/•OO 7 

smTTs / , o a 



= ^ / dz'z-'^^P°'^iz^,s) . (3.25) 



In the previous two expressions R means that we have to take the regular part of the large— z expansion of the 
integrand, whereas pole refers to the pole part at large z. 

The integrand functions have the asymptotic behavior H3.15|l . as shown in Appendix ^ we are, therefore, justified 
in applying the lemma discussed previously. Then, it is easy to see that 

Cr{s) = [/^]iog + si[I%,, I^m + 0{s) , (3.26) 

and 

Cp{s) = [P^'^'liog - s dz' In z^ ^' > + 0{s) . (3.27) 

The derivative can now be calculated easily and the previous results combine to give 

C(0) = [I%o, + [/P°'°]rcg - /P°'^(0) , (3.28) 

and 



cm = [/«],eg - /^(O) - / dzHnz^ ""' J'> 



(3.29) 







where we have anticipated the fact that [/P°''^]iog = 0, as shown in Appendix IbI The asymptotic expansion of /(z^, s) 
is also given in Appendix ^ where the various pieces appearing in (|3.28|) and H3.29|l are obtained. 

For C(0), [/^]iog, [/P°'1rcg and Ip°^%0) are needed, and Eqs. (jBTj and (jUl) provide them. Some algebra leads 

to the desired result: 

-1^0 ^'^(^+^--'^' . ("0) 

On the other hand, the evaluation of C'(0) requires the knowledge of [/^]iog, I^ {0) and of the integral piece appearing 
in Eq. (|3.29(l . These are calculated in Appendix IbI and the results are reported in formulas (|B9|I . I|B20|I and l)B38(l . 
which, combined together, lead to the following expression for C'(0)- 

C'(0) = C(0) ln(£sinh?7o)^ + Ai + A2 + A3 + A4 (3.31) 

where 

Ai ^ ^ - ^C'(0, 1/2) + ^C'(-l, 1/2) + ^C'(-2, 1/2) - ic'(-3, 1/2) + ^('(-4, 1/2) (3.32) 



137 d^(^ 

A2 = T 

7200 dx^ 



^ dxhix—-^+ f X ^(t>(x) dx , (3.33) 



120 Jo dx^ 
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A3 = - i / rfa; ^27ra: _ ^ ln[Z(ia; + 3/2)] - C.C. - - ln[Z(3/2)] , (3.34) 



73 251 T 31 1 

A4 = sinh^ Tin sinh^ Tin H cr^ sinh^ rjn H (t^ sinh^ r/n cr'* sinh^ rio , (3.35) 

1536 ' 3072 64 " 192 ' 48 ' ' ^ ^ 

where 

0(a:) = ^x-ia:3^1nZ(l/a:), (3.36) 

Z(x) = 2"r(l + a;)p-_"^/2(coshryo) . (3.37) 
Formulas H3.30|l and H3.31|) represent the main result of our paper. 

IV. LIMITING CASES 

The results obtained in the previous section for the zeta function and its derivative are valid for a scalar field of 
arbitrary mass m and coupling ^. Here we focus on the specific case of a massless conformally coupled field, as this 
allows us to compare our result to that of Refs. (2^i2^. By setting cr = 1/2 in Eqs. H3.30|l and ()3.31|) . the following 
expressions are found: 

1 17 

' ~ 6 1920 ' 

1 17 31 131 

C'(0) = ln(^sinhr7n)^ sinh rin sinh nn + Ai . 

^ ^ ' 6 1920 ^ ' ' 1536 ' 3072 



(4.1) 



These values can be compared with those of Ref. |2J| . After sorting out some transcription errors in Tables I and II of 
the mentioned reference, the result is found to disagree by a constant number. The question arises as to whether or not 
this difference is at all significant. Clearly, this constant difference can be reabsorbed by redefining the renormalization 
scale, and therefore it does not have any physical significance. However, from the mathematical point of view, the 
origin of such a difference is not so clear-cut. 

This disagreement has led us to consider another limiting case, which is obtained when the AdS curvature radius £ 
is large. This should reproduce the ball geometry, and corresponds, in our termin ology, t o £ >> tq, i.e. 770 << 1. This 
result has been computed by various authors using different techniques [sS, |^ 132113^ l33 | , and therefore it should 
provide quite a robust check on the result, and also an alternative derivation, although more involved than necessary, 
of the zeta determinant for the ball. Now, in the limit of t^q << Ij we find 



C(o) 



1 17 



6 1920 

1 17 , n 47 1 



C'(0) lnr2 + C'(0,l/2) (4.2) 

^ ^ ' 6 1920 ° 9216 64^ ^ ' ' ^ ^ ' 

+ lc'(-l, 1/2) + ^C'(-2, 1/2) - ic'(-3, 1/2) + ^C'(-4, 1/2) . 

This result is found in full agreement with those of Refs. [33. l33ll33 |. thus provi ding a robust check of our result. We 
note that such a hmiting case is not recovered neither by the result of Refs. [22,I3| nor by that of Refs. |2lll2^. 

In order to find some explanation for the difference, let us briefly reconsider the method used there. The first step of 
the method used in Refs. f2ll ||2^ I2M |2^ is a conformal transformation that changes the original background spacetime 
into a different one, where the evaluation of the zeta determinant is, in principle, easier. In Refs. |2ll 122 . 123ll24l | this 
corresponds to a coordinate transformation, dz = dr/a{r), upon which the original line element becomes 

ds^ = a'^{r){dz^ +dsl) . (4.3) 

In other words, a conformal rescaling of the metric by a(r) leaves us with a flat cylinder with a de Sitter cross-section. 
This is the starting point taken in the above mentioned articles. One immediately sees that the coordinate r lies in the 
range of < r < vq in the original frame, whereas in the conformally transformed one z lies between z{rQ) < z < 00, 
implying that the conformally transformed spacetime corresponds to a semi-inflnite cylinder and a more important 
point is that the conformal transformation is well defined at every point of the AdS bulk, except from the center where 
it breaks down. This observation raises the question as to whether or not this procedure is actually valid. Certainly 
it requires care. Obviously, in the two brane setup this problem does not exists since r is never 0. 
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V. 



CONCLUDING REMARKS 



The present article was devoted to providing an alternative derivation of the zeta function and of the functional 
determinant for a scalar field on a de Sitter brane background, which consists of a higher dimensional AdS bulk 
spacetime bounded by a de Sitter section. We considered the general case of a non-zero mass and coupling to the 
scalar curvature, thus generalizing previous results limited to zero mass and conformal coupling. 

For simplicity, we considered the case of a five dimensional bulk spacetime and Dirichlet boundary conditions. 
However, the result can be extended, with no additional technical problems and modulo some algebra, to other 
boundary conditions or higher dimensionalities. The choice of five dimensions was also motivated by the possible 
relevance of our calculation for the bulk inflaton model proposed in [ssl Is^ . 

One of the interesting points of the approach developed here lies in the fact that we do not make use of conformal 
transformations as is done in Refs. 122. 12 3. |24l and whose results, in the limit of large AdS radius docs not 
reproduce that of the ball given in Refs. [32, |33, yj] 

The basic tools of our calculation were a contour integral representation for the zeta function and asymptotics of 
the eigenfunctions, for which we have followed Refs. |27ll3lll3^ . In particular, the method devised in [23| has proven 
to be very powerful, and extremely useful in the case discussed here. 

As a check on the calculation, we have considered the limiting case of massless, conformally coupled fields, which 
was found to disagree with the result of [2^l23| . The difference amounts to a constant, which is physically harmless 
and can be removed by redefining the renormalization scale. This difference motivated us to consider another limiting 
case, which is obtained when the curvature radius of the higher dimensional AdS space is very large, leading to a 
ball-like geometry, for which extensive calculations of the zeta determinant are available In such a limit we 

recover those results and this should provide a robust check of our calculation, and, interestingly, this limit provides 
an alternative derivation, although technicall y u nnecess ary, of the ball determinant. We note that the results for the 
derivative of the zeta function given in Refs. |2ll 123. 12M l24| do not recover this limit. 

Various generalizations of the work presented here are possible. Extensions to higher dimensionalities, different 
boundary conditions, two-brane setups and higher spin fields should follow without additional difficulties and only 
a larger amount of algebra might be needed. Looking at the possible relevance of these kinds of calculations in the 
bulk inflaton model and more generically in brane world cosmology also deserves further study. We hope to report 
on these issues in our future work. 
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APPENDIX A: UNIFORM ASYMPTOTIC EXPANSION OF THE LEGENDRE FUNCTIONS 

Following Refs. [s^l and [s^, the uniform asymptotic expansion of P~^^^^{cosh.ri) can be obtained. In this appendix, 
we simply outline the procedure and extend the results to the case of the logarithm of the Legendre function (cosh 77) . 

For large values of /x, the solution of the Legendre differential equation can be written in the well-known WKB form. 
When substituted into the original equation, this leads to a set of recursive equations that allow one to determine the 
expansion. Omitting the details, the result is 
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where the functions t and ^ are given by 



1 



(A2) 



t = 



\/l + sinh^ 77 



and 



^ = In 



T sinh T] 



T tanh ^ tanh 

T 



1 



1 



+ 1 . 



(A3) 



■\/l + sinh^ 7] + cosh 77 



Tt cosh T] 
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We have defined t = n/v and the coefficients ai can be computed recursively. They are quite lengthy for increasing 
values of i and since they are not used directly we do not report them. The interested reader is addressed to the work 
of Refs. 1^3, Ell, where they are also derived. 

It is a relatively simple task to get the logarithm of the previous expression and this can be achieved by using a 
symbolic manipulation program. The result is 



InP 



'^_^i/2(cosh77) - In 



n=l 



where it is essential to note that the coefficients En{t, t) are bounded in the full range Q < t < oo and for large t they 
scale as inverse powers of r. For — > oo they exhibit a power law growth of finite order in v. This can be checked for 
the first four coefficients i?„, which are found to be 

El = -4 (Sr^i^ cosh^ 77 - 3i cosh 77 + 1 - Sr^t cosh 77) , 

24(t^ ^1) 

E2 = — r4- ^ {-l + t^ cosh^ T] + Tr^t^ cosh^ 77 - ~ QtH^ cosh'' 77 + r^t^ cosh^ 77 - Gr^i"* cosh* 77 

i6(t^ — ly ^ 

+?>tH^ cosh^ 7;) , 

E3 = ^ -T (ll888lT*t^ cosh^ 7? - 51765x2^3 cosh^ 77 + 18270x2^ cosh 7/ - hnmrH^ cosh^ 77 

40320(r^ — l)'^ ^ 

+2835x4* cosh 7/ + 33453x4^ cosh^ 77 - 2625x^*3 cosh^ 77 - QmihrH'' cosh^ 77 + 33453x^*5 cosh^ 77 
-mQlhrh'^ cosh^ 77 + 38675r'^*^ cosh^ 77 - 98 + 2835i cosh 77 - 294r2 - 2625*^ cosh^ 77) , 

Ei = r4 7T (565r®ii2 ^^^^12 ^ isQSr^i^ cosh* 77 - 2282T^i^ cosh^ 77 - 2%StH^ cosh^ 77 

128(r2 - 1)4 ^ 

-1356r^t*° cosh*" 77 - 288T2i^ cosh^ 77 + 542T^t* cosh* 7; - 2282r*t^ cosh^ 7; + 542r^t* cosh* 77 
-18r^t^ cosh^ 77 - 284T^t^ cosh^ 77 - 284r**t^ cosh^ 77 + 13T*i* cosh* 7; + 1062r^t** cosh® 77 
-1356T®t*° cosh*" 77 + 5 + 30r2 + 5t* + 13t* cosh* 77 - lU"^ cosh^ 77 + 1062T*i® cosh® 77 
+3114T^t®cosh®77) . (A5) 

By inspecting the previous expansion and recalling the above-mentioned properties of the coefficients Ek (t, t) , one 
notices that it has the same structure as Eq. (|3.15|) . 



APPENDIX B: EVALUATION OF I^{z^) AND 1"^°^" {z'^) . 
1. Evaluation of [/^]iog- 

From Eq. dsn, we can find the various pieces that appear in Eqs. H3.28|l and (|3.29(l . Let us first consider [/^]iog. 
From Eq. ljA4|l . we see that the coefficient of the logarithmic piece, for large z, is 

To deal with this sum, we introduce a generalized C-function: 

C.(s)-Ed,i.-« (B2) 

j 

which can be expressed, in five dimensions, in terms of Hurwitz ^-functions. Trivial manipulations give 

[Aog = -^C.(2s-l)-ic.(2s) , (B3) 

which, in five dimensions, becomes 



11 



2. Evaluation of IP°^''{z^). 



The term /p°''^(z^) is more involved to evaluate. We need to rewrite the uniform asymptotic expansion ljA4|l in 
terms of inverse powers of v and then perform the summation over j. This will allow us to express, in five dimensions, 
expansion ljA4p in terms of Hurwitz ^-functions from which the pole part can be extracted. The calculation is rather 
lengthy, although straightforward. Here, we simply quote the result, which can be written as follows: 



where we have defined the following quantities for notational convenience: 

sinh^ ryo 



UJi{x) 

UJ2ix) 
UJ3{x) 

UJ4{x) 



16(1 + .x^ sinh rjo) 
sinh'^ 770 
16(1 + x'^ sinh'^ 770)'^ 
1 



(4 + (7 - lOx^) sinh^ 770 + x^{-8 + x^) sinh^ tjo) , 
(1 - 42:^ + a;^(-4 + x^) sinh^Tyo) , 



[ sinh^ 770 ( 8 - 32x2 + sinh^ 770 (13 - 2Ux'^ + 288x"* 



128(1 + a:2 sinh^ 770)6 
- 2a;2(116 - SlSx^ + 1162:'^) sinh^ 770 + a;^(288 - 244a;2 + 13a;'') sinh^ 770 
+ 8a;6(-4 + a;2)sinh%o ) ] , 

sinh^ 770 / 1 cr^ sinh^ 7/0 



1 + x'^ sinh 770 V 2 1 + x'^ sinh 770 



(B5) 



(B6) 



The absence of logarithmic terms in implies that [/P°''']iog ~ 0. From the expression ljB5|l . one easily finds that 



j-jpolcj 



rog 



lim /P°"=(z2) = 



(B7) 



and 

P-^m = (1 + 2smh2 ,0) - 1/4) + ^ (.2 - 1/4)^ 

both of which are required for the evaluation of C(0)- 

The integral of the pole piece is readily evaluated and the result we find is 



dz In ; 



,d/P°l'=(z2) 



47 



9216 1536 



73 . ,2 251 . . 

smh rjo + smh 770 



3072 



1 



(B8) 



7 2 • 1 2 2 • 1 4 ^ 4-14 

— a smh 77n <J smh rin a smh 77n 

64 '° 192 48 

sinh' 770(4^2 - l)(-6 + sinh' 77o(-13 + Aa^)) ln(£ sinh 770) 

(Do 



(B9) 



3. Evaluation of /" (O). 

To evaluate /^(O), we follow Ref. [27| once again and employ a more expedient approach based on the Abel-Plana 
summation formula. For z = we have 

^(0,^)-E4^' (BIO) 

3 J 

fiiy,) = i _ 1^^.^ 2''^r(l + 7.,) P;:i (cosh 770)] . (Bll) 

The validity of the form of the function, fii'j), which we use in the Abel-Plana summation formula is discussed in 
[27|. By applying the Abel-Plana formula, which allows one to convert the sum (|B10|) into an integral, we get 



1(0, s)= -^—dx + i — dx + -f{iyo), (B12) 
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where we retain the regularizing factor, Vj^, only in the first term, because all other terms are finite as s ^ 0. The 
only non-trivial term to compute in the expression is the first one. In order to evaluate it, we split the integral into 
three pieces, for convenience. Some algebraic manipulations give 



^.^^ 



dx 



x^'-''(f){x)dx + / x-'=<j){x)dx- / 



2/3 



where we have put s = where possible. The function 0(x) is given by 

2i/"r(l + l/a;)P"_'("(cosh77o) 



{x) = ^(x~^x^\ In 



(B13) 



(B14) 



Integrating the first two pieces by parts six times, we obtain 







dx — 



X 



2s-5 



J2s-i 



2s- 5 
1 



^{x) 



d(j){x) 



,2s 



(2s-5)(2s-4) dx 



(2s-5)(2s-4)---2s 



(2s-5)(2s-4)---2s dx^ 
1 



dx 



X ^ X ^ d(f){x) 



In a; d^(f){x) 
-120 dx5 



120 7i ^"^^^^'^^"X/a" '^^''^ 



1 

(B15) 



This can be appropriately analytically continued to s = 0, giving 



f{v,) 1 , 137 



Summarizing, we find 



where 



240s dx^ 



7200 



In a; — f-^^dx 
dx^ 



2/3 



. (B16) 



rZ%)d^^i/poie(o)+/R(o) 



/P°''=(0) 



1 d^(j) 
240 



(B17) 



(B18) 



x=0 



/?(0) 



137 d^c?!) 



7200 dx^ 



1 

/ dxlnx — - 

-_o 120 Jo dx^ 



X ^<j)[x) dx 



2/3 



(B19) 



A nice consistency check on the previous evaluation is given by the fact that the expression for /p°'°(0), 
evaluated previously using the asymptotic expansion agrees with (jB18|l . 
The previous results can be combined together to get 7^ (0). We find 



/H(0) = ^ 
^ ' 7200 dx^ 



-^-^ I dxlnx-j^ — I X ^0(x) dx 



.0 120 7o 



2/3 



+ i I dx 
'o 



(ix + 3/2)3/3 - («x + 3/2)/12 



+ iln[23/2r(5/2)F-_\//Vcosh,yo)] ■ 



1/2'" 



ln[2"+3/2r(zx + 5/2)P^"_:^^2'^'(coshr;o)] + c.c. 



(B20) 



4. Evaluation of [/"] 



rcg- 



We make the final effort to obtain the regular part of /^(z^). From previous arguments, we understand that such 
a piece comes from the terms in the asymptotic expansion that scale as z'^ , which then we have to sum over j. Thus, 
from I A4l we have, for s ^ 0, 

[/«]reg = - J(ln27r + ln(£sinh77o)) djilnu^ + lnT{iyj)) , (B21) 
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where we have used the fact that 



Urn Vd,!/]-^'^ = 



(B22) 



The first two terms can be computed easily, whereas to deal with the last sum in ljB2ip we proceed as follows. First 
we use an integral representation for the logarithm of the F— function [33l Isft] . which gives 

(B23) 



dj (In i., + In F(i.,)) = i In 27r ^ + A + 



where we have defined 



(B24) 



and 



1 1 



2 t e* - 1 



. dt 

T 



(B25) 



Let us first deal with 9. ft is possible to sum the series appearing in this expression by using the relation 

-t/2 



e* - 1 



(B26) 



Differentiating this relation one and three times, one immediately arrives at 



1 



1 1 



3 7o V2 t l-e-\ 
which can be expanded to write it in terms of elementary integrals: 



1 / d 

4 \dil- e-* 



dt 
T 



^£-2g~9t/2 ^e-2g-7t/2 ^e-2g-5t/2 ^e-2g-3t/2 

"^(l-e-*)4 ~ ^(l-e-*)3 ~ '*(l-e-*)2 " (1 - e"*) 

-llt/2 ie-le-9t/2 13i.~lg-7t/2 ^^-Ig-St/Z 1^6-1-34/2 

'2— ^+6- 7TT + -rTZ 7T^+3' 



(l-e-*)5 (l-e-*)4 2 (l-e-*)3 (1 - e-*)2 2 (1 - 



(B27) 



(B28) 



Here we have introduced a regulating factor e, which ensures the convergence of the expression for 3?e > 2. The limit 
e ^ will be taken at the end. 

All the above integrals can be evaluated starting from the standard formula 



<^i(a,6,/3) 



dt 



1 



-F(l + a)C(l + a,V/3) , 



/o 1 - e-P* t-" ^1+° 
which, by repeated differentiation with respect to /3, produces the following relations 

^oo g-(&+/3)t g 



(B29) 



(B30) 



ip3{a,b,f3) 



ip4ia,b,P) 



(p5{a,b,P) 



oo g- (6+2/3)4 



(1 




-2-a 


e" 


-(6+3/3)t 




(1 






e 


-(f)+4/3)t 





1 Id 

--^2(a + 1,^/5) - 2a^'^2(a,6,/3) , 

2 19 



(1 - e-/5*)5 <-4-a 4 



1 9 

934(0 + 1, b, (3) - |^V4(a, b, 13) 



(B31) 



(B32) 



(B33) 
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The previous relation (|B28|I can then be expressed in terms of the functions tpj (a, h) = (/Sj (a, 6, 1) . Simple calculations 
lead to 

e = -2994(e-5,3/2)-5^3(e-4,3/2)-4(/j2(e-3,3/2)-^i(e-2,3/2) 

+2(p5(e - 5, 3/2) + 6¥>4(e - 4, 3/2) + ^-^^{e - 3, 3/2) + 3¥>2(e - 2, 3/2) + ]^ifi{e ~ 1, 3/2) . (B34) 
The limit e — > can now be taken and the result is found to be 

^ = - ^C'(-2, 1/2) + ^C'(-4, 1/2) . (B35) 

The term A can be evaluated, by using the following relation: 

"^■E-r = -E-ri--.- (^36) 



ds 

J 



A straightforward computation leads to 



A = -^C'(-3, 1/2) + ^C'(-l, 1/2) - ic'(-4, 1/2) + ^C'(-2, 1/2) . (B37) 
Some simple algebra allows us to combine the previous results to arrive at 

[^""^^^^ ^ Y 9^ ln(£sinh,7o) - ^C'(0, 1/2) + ^C'(-l, 1/2) + ^C'(-2, 1/2) - ^('(-3, 1/2) + lc'(-4, 1/2) . (B38) 
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